An analysis of slippage effects on a solid sphere enclosed by a non-concentric cavity filled with a couple stress fluids

This investigation shows the effect of slippage on the slow spinning of a rigid sphere covered by a non-concentric spherical hollow full of an incompressible couple stress fluid. Moreover, the velocity slip conditions are employed on surfaces of both the rigid sphere and the cavity. In addition, the solid sphere and the cavity are rotating axially at various angular speeds. The solution is obtained semi-analytically at low Reynolds numbers utilizing the superposition with the numerical collocation approach. This paper discusses the hydrodynamic couple exerted by the fluid on the internal particle. The dimensionless torque increases as the slip and spin slip increase by 99%, the couple stress parameter by 49%, and the separation parameter by 79%. Additionally, the non-dimensional torque decreases with the increase of the size ratio by 89%. Consequently, it is found that all the results agreed with the corresponding numerical analysis in the traditional viscous liquids and the revolving of two eccentric rigid spheres with no slippage (Al-Hanaya et al. in J. Appl Mech Tech Phys 63(5):1–9, 2022).


Mathematical formulation
By the notion of a low Reynolds number, without the presence of body forces and body couples, the field constraints controlling the steady motion of an incompressible couple stress liquid are dictated by 2 : Here, the constant µ is the fluid viscosity, η is the viscosity of 1st couple stress, u is the velocity vector of the fluid, and p is the fluid pressure.If the relation (2) tends to the classical equation of Navier-Stokes.The tensor forms of t ij and m ij are 2 where η ′ is 2nd couple stress, m is a tensor trace of couple stress, the two tenors are Kronecker delta, e ijk is alternating tensor, d ij is deformation rate tensor, and ω is vorticity vector, the two last concepts are formed as: The enforced boundary constraints may be used to derive the scalar quantity that was mentioned in relation (4).Additionally, one may explicitly define it as 1 by using the second relation of ( 4) and the concept of (5): The physical constants in the fundamental Eqs. ( 3) and (4), as well as the equation for movement (2), are presumed to adhere to the following constraints in 2 : (1) (3) Furthermore, suppose that the surface of the sphere is subject to the following circumstances., r = a: a. Slippage restriction where is the slippage parameter changing its values from zero to infinity.This coefficient is only related to the type of fluid and the material's surface.Furthermore, the perfect slip situation becomes possible when the slip coefficient disappears, and the traditional no-slip case may be inferred as a specific instance in this study when the slip parameter approaches infinity.The slippage condition of the boundary has recently been used to solve several viscous fluids 9-12 and micropolar fluids 26,27 issues.b.The prevailing condition is the absence of couple stresses 2   where is the unit normal to the surface of the solid sphere.Stokes 1 has suggested the boundary conditioning Eq. (9).Only in this situation can mechanical interactions at the borders produce a force distribution., according to physical theory.

Solution of the problem
Assume that the rotational movement of a spherical object of radius moves symmetry about its axis within an incompressible couple stress liquid.The spherical systemic procedure is established at the center of the sphere, the field functions are not dependent on φ .Further, the velocity and vorticity vectors are represented by: Substitute Eq. ( 10) into the momentum Eq. ( 2) by eliminating the pressure, the subsequent p.d.e is obtained as: where the material constant 1/κ = η/a 2 µ , is taken into consideration as a polarity indicator for the couple stress fluids approach, and the Stokesian indicator of axial motion is: Moreover, from Eq. ( 5) the non-vanishing vorticity components ω r and ω θ , are: Furthermore, the tangential stress is calculated by El-Sapa and Almoneef 7 as: We obtain the following couple stresses by using the tensor relation ( 4): www.nature.com/scientificreports/ The differential Eq. ( 12) has the following generalized solution: where the two functions and K n (.) the first and second types of modified Bessel functions of order n , respectively.Also, P 1 n (.) denotes the corresponding Legendre polynomials of order the first type.Applying Eq. ( 18) to Eqs. ( 14) and (15), the vorticity components are obtained as: The couple is determined by applying Eqs. ( 19)- (20) to Eq. ( 17): The boundary conditions ( 9) can be written as (17e) (κr) Employing the obtained Eqs. ( 18), ( 21)- (25), and ( 27) into ( 16), we get: where As a result of the axisymmetric particle being impacted by the fluid flow, a torque is generated that has a magnitude of 28   Hydrodynamic interaction of a rigid sphere enveloped by a spherical cavity filled with a couple stress fluid For this simulation, it is assumed the annulus between the solid sphere, a 1 and a spherical cavity, a 2 is filled with a constant density of couple stress liquid.Therefore, the sphere and the spherical cavity are rotating around a connecting line of its centers with distinct angular speeds 1 ad 2 , respectively and at a distance h from their centers as shown in Fig. 1.Consider that u θ are the components of velocity and vorticity as a result of the presence of the solid particle a 1 without the spherical cavity a 2 and u θ are the components of velocity and vorticity of the spherical cavity a 2 without the solid particle a 1 as shown in Fig. 1.Additionally, the subsequent relations link the coordinate systems (r 1 , θ 1 ) and (r 2 , θ 2 ) together as: Hence, the boundary conditions are linear so the principle of superposition can be applied.Thus, the field functions are represented as: (26) (κr) (κr) (κr) Vol:.( 1234567890) From (32) and use ( 18), ( 22), (28) the field functions q φ , ω r , ω θ , m rθ , t rφ are: www.nature.com/scientificreports/Accordingly, by applying the boundary conditions ( 34) and (35) to Eqs. ( 35), (38), and (39), we get the following system: (38) (39) (κr 1 ) (κr 1 ) (κa 1 ) where These constants A n , B n , C n , D n obtained by solving 4N simultaneous linear algebraic Eqs.(40)-(43) provided by the infinite series has been truncated to N terms.To satisfy the boundary criteria at a limited number of discrete locations on the generating arcs of the spherical boundaries, the boundary collocation technique will then be used.The desired unknowns A n , B n , C n , D n are then determined by numerically solving the resultant system of equations using the Gauss elimination technique.On the semi-circular longitudinal arc of each particle surface from θ = 0 to θ = π , the collocation technique (Ganatos et al. 1980) applies the boundary conditions at a finite number of individual points and reduces the infinite series in Eqs. ( 40)-(43).The coefficients matrix becomes unique if these points are employed, as shown by looking at the system of linear algebraic equations for the unknown constants A n , B n , C n , D n .The strategy suggested in the literature, such as that used by Ganatos et al. in 1980 to choose the collocation points, is what we apply to avoid this singular matrix and obtain high accuracy: Four fundamental collocation points on each spherical particle are taken at θ i = ε, π/2−ε, π/2+ε, π − ε on the half unit circle 0 ≤ θ i ≤ π at n any meridian plane, where ε is provided by a minimal number to prevent the singularity at θ i = 0, π/2, π .The other points are chosen as mirror-image pairs with θ i = π/2 and are uniformly spaced around the two-quarter circles, omitting those singularities.The linear algebraic equations are solved (41) (ℓa 1 ) (ℓr 1 ) (ℓr 1 ) D n P 1 n (ζ 1 ) its axis.It is also known as the moment of force.Physically, speed is defined as the amount of distance traveled in each amount of time, but angular velocity refers to the rate of rotation of the body and the number of revolutions in each amount of time.Thus, the non-dimensional torque is inversely proportional to the angular velocity.Consequently, that the torque decreases with increasing, ω which agrees with the physical concepts.in addition, the torque slowly decreases for the positive values of the angular velocity with the increase of size ratio but for the negative values, it changes its direction to up. Figure 2b displays the torque for various values of slippage parameters for fixed the cavity and the solid sphere rotates with value one.Hence, the torque increases with the increase of both the slippage parameters and the size ratio.This mode is like Motor Mode, the induction motor torque swings in this mode of operation as the slip changes, going from zero to full load torque.From zero to one is the slide.At no load, the value is 0; at rest, it is 1.But for the fluid, the slip parameter varies from zero to infinity, where zero denotes perfect slip, the values in between are the partial slippage, and infinity denotes the no-slip condition, the last is the limiting situation for the work of Amal et al. 6 .The curves show a clear relationship between the torque and the slip.In other words, the amount of torque produced increases with slippage and vice versa.Additionally, Fig. 2c differs from Fig. 2b with the value of the separation parameter, and also the cavity has a partial slip which makes the torque diminish rapidly with the increase of the size ratio, the torque value appears to be minimal when the particle is in a concentric position inside the cavity ( δ ≈ 0 ) as expected.
We have found that our placement results of the torque in a concentric position are very similar to the analytical solution available in the literature.Moreover.Figure 2d for no-slippage indicated the advancement of torque with the improvement of the first couple stress parameter.Furthermore, Fig. 3a-d illustrate the distribution of torque against the velocity slippage on the solid sphere for certain values of the rest parameters.Figure 3a displays the torque growths with the growth of the slip parameter on the cavity and at the same time increases with the slippage on the solid sphere.On the other hand, in Fig. 3b the torque diminishes with the increase of the size ratio, a 1 /a 2 and the torque as mentioned previously increases with the slip parameter.Figure 3c  For various values of the indicated parameters, the distribution of torque on the sphere versus the size ratio with (a) η = 0.01, η′ = 0, a 1 /a 2 = 0.5, δ = 0.01, ω = 0 , (b) η = 0.01, η′ = 0, β2 = 0.5, δ = 0.01, ω = 0 (c) η = 0.01, η′ = 0, β2 = 0.5, δ = 0.01, a 1 /a 2 = 0.5 , (d) ω = 0.5, η′ = 0.5, β2 = 1, δ = 0.01, a 1 /a 2 = 0.5.
of slippage parameters and the torque declines with the increase of the angular velocity ratio as usual.Thus, in Fig. 3d exposits the torque has more significant for the improvement of the first couple stress parameter and agree with the limiting case of viscous fluids.Table 3 shows the exact value was taken when N = 90 for the convergence of the normalized torque with various parameters.
Finally, the distribution of torque versus velocity slippage on the solid sphere for fixed values of the pertinent parameters is presented in Fig. 4a and b as a result of this investigation, the torque increase with the increase of the second couple stress parameter increases.Additionally, Fig. 4b and c expressed the torque versus the separation parameter where the improvement of torque has been shown in Fig. 3b by increasing the distance from the centers between the solid sphere and the spherical cavity.Therefore, in Fig. 4c the torque is directly proportional to the slippage parameter and this relation affects the torque with the separation distance that the torque inclined with the growth of the separation parameter.

Conclusion
In this research, we study the interfacial slippage effect and the steady incompressible rotation of a couple stress fluids around a rotating sphere.Therefore, the graphs are used to give a numerical analysis of the torque operating on the inner solid sphere's surface.As a result, raising the couple stress coefficient results in an expected increase in torque.Additionally, it is established that the torque is significantly influenced by the second viscosity parameter.It elevates the torque's worth.In addition, it is shown that the size ratio and separation parameter have more significant on the couple stress fluid flow, especially for small values.Finally, it is determined that the slip parameter has a significant influence in raising the torque value.The motivation for studying flow slip boundary conditions comes from the possible applications in a variety of engineering and applied scientific fields, as well as from a serious grasp of hydrodynamics, which serves as the theoretical basis for the design and construction of nanofluidic devices.Additionally, the development of shale reservoirs depends heavily on a knowledge of slip flow behavior in the nano-porous medium.The future study can be applied to this work in the effect of permeability  of porous medium and magnetic field such as in 33 and 34 .Additionally, the impact of oscillation, the fractional approach and electro-osmotic can be employed in this study such as 35-37 .
Table 1.The non-dimensional torque exerted on the internal sphere in this situation when the exterior sphere is stable and the interior sphere is spinning, with zero second couple stress and η = 0.01.Table 2.The non-dimensional torque exerted on the internal sphere in this situation when the exterior sphere is stable and the interior sphere is spinning, with zero second couple stress. T

Figure 1 .
Figure 1.The geometrical shape of a sphere covered by a cavity filled with couple stress fluid.